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Abstract: 

In the semi-classical treatment, i.e. in a classical black hole geometry, Hawking quanta 
emerge from trans-Planckian configurations because of scale invariance. There is indeed no 
scale to stop the blue-shifting effect encountered in the backward propagation towards the event 
horizon. On the contrary, when taking into account the gravitational interactions neglected in 
the semi-classical treatment, a UV scale stopping the blue-shift could be dynamically engendered. 
To show that this is the case, we use a non-perturbative treatment based on the large- N limit, 
where N is the number of matter fields. In this limit, the semi-classical treatment is the leading 
contribution. Non-linear gravitational effects appear in the next orders and in the first of these, 
the effects are governed by the two-point correlation function of the energy-momentum tensor 
evaluated in the vacuum. Taking this correlator into account, backward propagated modes are 
dissipated at a distance from the horizon tx Gn when measured in a freely falling frame. (G is 
Newton's constant and k the surface gravity.) This result can be also obtained by considering 
light propagation in a stochastic ensemble of metrics whose fluctuations are determined by the 
above correlator. 



*An earlier version of this work was published in Proceedings [21] in 2002. We postponed submitting it to 
the arXiv in the hope of improving the evaluation of radiative corrections. We have recently modified the text, 
corrected some mistakes, and added remarks on this difficult point which still needs work - that hopefully someone 
will take over. 

^ E-mail: Renaud.Parentani@th.u-psud.fr 
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1 Introduction 



In his original derivation [T] , Hawking considered the propagation of quantum radiation in a fixed 
background metric, that of a collapsing star. This means that the metric is once for all deter- 
mined by the energy of the collapsing star. It is therefore unaffected by the quantum processes 
under examination. In this approximation, the radiation field satisfies a linear equation (in the 
absence of matter interactions). One then finds that in- falling and outgoing field configurations 
are completely uncorrelated near the black hole horizon^] This is explicitized by the fact that 
the connected part of the two-point correlation function of T uu and T vv , the energy-momentum 
tensor of outgoing and infalling configurations, 

{T vv (x) T uu (x')) c = (T vv (x) T uu (x')) - {T vv (x)) {T uu (x')) , (1) 

vanishes in the (Unruh) vacuum. Nevertheless, Hawking radiation is pair creation. This is 
perfectly consistent with eq. ([T]) since the pairs are composed of two outgoing quanta, one of 
each side of the event horizon. The external ones form the asymptotic flux whereas their partners 
propagate towards the singularity at r = 0. Upon tracing over these inner configurations, one 
gets an incoherent flux described by a thermal density matrix[3jH 

From this fixed background description, one may go one step further by performing a mean 
field approximation, i.e. by including the metric change determined by Einstein's equations 
driven by the expectation value {T^ u ). One then finds that this expectation value is regular 
[TJ [U [3]. This guarantees that the black hole will adiabatically evaporate while keeping 
the regularity of the near horizon geometry. This regularity in turn implies that the infalling 
and outgoing configurations will stay uncorrelated. Therefore, the correlation function ([1]) still 
vanishes in the semi-classical treatment. 

This adiabatic evolution would provide a reliable starting point for including perturbatively 
radiative corrections were another feature of black hole physics not present. Namely, the field 
configurations giving rise to Hawking quanta possess arbitrary high (trans-Planckian) frequencies 
near the horizon [9jllO } [lT j l3"]. When measured by free falling (FF) observers at r, the frequency 
of an outgoing quantum of asymptotic energy A grows as 

ft K ~TTT~ ■ (3) 

This implies that wave packets centered along the null outgoing geodesic u had FF frequency 
growing as f2 oc \e KU when they emerged from the collapsing star. Unlike processes characterized 



In this paper we consider only s-waves and neglect the residual potential barrier which partially back-scatter 
some Hawking quanta thereby inducing some correlations. These however play no role in what follows and shall be 
ignored. In the next equations, v and u are radial advanced and retarded null coordinates. In the Schwarzschild 
geometry, they are given by v = t + r* , u — t — r* where r* = r + 2M ln(r/2A/ — 1) is the tortoise coordinate [2]. 

2 The non- vanishing character of the correlation function (T uu T uu )c across the horizon demonstrates that the 
partners of Hawking quanta are inside outgoing configurations. An explicite calculation gives[U[S] 

(T„„(w) T uu (u' L ))c oc \u — u'l + iTv/n\~ 4 , (2) 

where kul = ]xi(kUk) is a null coordinate for the inside configurations. Uk is the usual Kruskal retarded time; 
it vanishes on the horizon, and is positive inside the hole, k = 1/4M is the surface gravity. (It fixes Hawking 
temperature Th = k/2tt. We work in Planck units: c = fi = Mpi anc k = 1.) The smooth maximum of this 
two-point function for opposite points, i.e. u = u' L or Uk = —U' K , is a direct consequence of the fact that each 
pair is indeed composed of two outgoing quanta leaving on either side. 
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by a typical energy scale, the relation oc Xe KU shows that black hole evaporation rests, in this 
scenario, on arbitrary high frequencies. This conclusion drawn from the analysis of wave packets 
is confirmed by the study of the non-diagonal matrix elements of T au . These characterize the 
fluctuations of the flux around its mean value. As shown in [U [TTj [3], unlike the expectation 
value (the diagonal part) which is regular and of the order of M -4 , these matrix elements are 
singular on the horizon, i.e., their Fourier content is characterized by FF frequencies which grow 
according to eq. ([3|). 

This growth is not an artifact of a bad coordinate choice. Indeed, as emphasized by 
't Hooft [12j, and explained below in Section 3, gravitational interactions between the con- 
figurations giving rise to Hawking quanta and in-falling quanta also grow according to eq. (|3|). 
This questions the validity of the semi-classical treatment and the vanishing of eq. ([1]) H 

In questioning the validity of the semi-classical description, two issues should be distin- 
guished, see e.g. section 3.7 in [3J. First, there is the question of low frequency O(k) modifica- 
tions of the asymptotic properties of Hawking radiation, and secondly, that of high frequency 
modifications of the near horizon physics. Since all thermo-dynamical reasonings indicate that 
the asymptotic properties (namely thermality governed by k and stationarity) should be pre- 
served, the problem is to reconcile this robustness with the radical change of the near horizon 
physics which is needed to tame the growth of gravitational interactions. This is not an easy 
problem: Indeed, the perturbative analysis of near horizon interactions performed in [l3j leads 
to back-reaction effects which grow like O in eq. ([3]). This threatens the stationarity of the 
flux and therefore questions the choice of the treatment which is adequate to go beyond the 
semi-classical approximation. 

As a first step towards a full quantum gravitational treatment, inspired by [E]-[TH], we 
proposed[19l I20| I21j a non-perturbative treatment of the interactions occuring in the FF vacuum 
based on {T vv T vv )c, the correlation function of infalling configurations. As discussed in more 
detail in what follows, this treatment emerges in a large N limit, where N is the number of 
copies of the quantum field. In physical terms, in this limit, in-falling configurations act as an 
environment for the outgoing quanta and their gravitational interactions express themselves in 
terms of a stochastic ensemble of metric fluctuations. The statistical properties of the latter 
are determined by the fluctuations of the N fields in the FF vacuum. Their main effect is to 
dissipate the outgoing modes traced backwards towards the horizon at the locus where their FF 
frequency ([3]) reaches 1/(7, i.e., at r = 2M + a for modes of Killing frequency equal to k. 

In our model, the gravitationally induced length-scale0 is 



rr\T^ y/N]n(A/K). (4) 



M 2 

m Planck 



In the above A is a high frequency cut-off whose value requires further study to be fixed. It is 



3 In section 9 of the review [12] . one reads "Any decomposition of Hilbert space in terms of mutually non- 
interacting field quanta will be hopelessly inadequate in this (near horizon) region." 

4 Added comment. When presented in meetings, this result was skeptically received on the basis that a is much 
shorter than the Planck length for macroscopic black holes. Indeed one has a oc L% lanck /GMi,h- However it should 
be pointed out that a, which is the proper distance between 2M and 2M + a measured in a FF frame, actually 
corresponds to a proper distance of the order of Lpi anc k when measured along surfaces of constant (Killing) time. 
It should be also recalled that r is defined by the square root of the area of 2- spheres/ An , and is therefore coordinate 
invariant. On the contrary, the notion of distance from the horizon depends on the choice of the time slices used. 
In this sense our effective propagation law does not singles out a prefered frame. This should be compared with the 
procedure involved when using an a priori given dispersion relation, see \22 )j - \27\[ . I am grateful to Ted Jacobson 
for discussions on this. 
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important to notice that in spite of these dissipative effects, the gravitational interactions do 
not significantly affect the asymptotic properties of Hawking radiation. As shown in [TB], the 
asymptotic corrections scale indeed as (net) 2 , hence they are negligible for macroscopic black 
holes. It should also be noticed that the dissipative effects break the 2D local [9] Lorentz 
invar ianceH 

An unsolved question concerns the range of validity of our treatment. This is a complicated 
question whose final answer requires a better understanding of quantum gravity. Let us never- 
theless make some remarks. First, this question closely follows that concerning the validity of 
the semi-classical treatment Secondly, our analysis indicates that the semi-classical treatment 
fails before our treatment. 'Before' should be understood radially, given the blue shift effect 
([3]) encountered in the backward propagation of outgoing configurations. What emerges is a 
kind of Russian doll structure in which quantum gravity progressively dominates the physics. 
Far away from the hole (r — 2M 3> 2M) one has outgoing thermal (on shell) radiation. In 
a first intermediate regime (a <C r — 2M <C 2M) the propagation of outgoing modes is still 
governed by the d'Alembertian but observers at fixed r and free falling ones perceive quanta 
differently (in that the locally defined FF vacuum no longer agrees with the Killing (Boulware) 
vacuum). It is in this regime (well described by the semi-classical treatment) that Hawking 
radiation gets established, see eqs. (76-83) in [5]. This description based on modes stops to be 
valid when reaching a Jacobson's time-like boundary jlOj . at r ~ 2M + cr, when outgoing modes 
get progressively entangled to the infalling configurations, thereby loosing their 'mode' quality. 
The principle aim of this paper is to analyze this transitory regime. Deeper in r, one has some 
unknown regime governed by Planckian physics. This physics presumably also occurs around 
us but stays well hidden inside its Planckian husk in the absence of a good microscope. 

Below we discuss two aspects. We first explain why we adopted a treatment based on the 
large N limit, and then we further discuss the relation between this treatment and the dispersive 
physics involved in acoustic black holes. These two subsections can be skipped in a first reading. 

1.1 The choice of the treatment: the large-A*" limit. 

In order to compute/estimate quantum gravitational corrections to Hawking radiation in the 
absence of a theory of quantum gravity, one should adopt an approximative treatment allowing 
to compute radiative corrections to some physical quantities. In this paper, we have chosen a 
statistical treatment based on a large- N limit, where ./V is the number of copies of the quantum 
matter field. The reason for this choice are the following. 

First, the semi-classical treatment is the leading contribution in the large- ./V limit. This 
is most simply understood in a path integral approach [29j: By performing a saddle point 
approximation in the evaluation the one-loop effective action, one verifies that the location of 
the saddle is determined by the semi-classical Einstein equations. Hence, the semi-classical 

5 In the vicinity of a black hole horizon, free propagation of s-waves is governed by a 2D Lorentz (and scale) 
invariance in the u, v plane, see Section 2. When including radiative corrections, the dressed propagator looses 
this property. This is similar to the fact that the self-energy of an electron immersed in a thermal bath of photons 
is not Lorentz invariant either 28 . In both cases, the integrands governing loop corrections are not Lorentz 
invariant, more on this in Section 5. 

6 The validity of the semi-classical treatment has been often questioned in rather general terms. However, a 
significant answer requires to find physical quantities (i.e., matrix elements of observables) which are incorrectly 
evaluated in this treatment and to propose improved expressions for the same quantities in order to see the 
discrepancy. We shall provide an explicit exemple in Section 4. 
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treatment is the mean field approximation (Hartree) in which the N copies of the radiation field 
propagate in a self-consistent classical geometry: a solution of Einstein equations driven by N 
times (T^), the mean value of the energy- momentum tensor of one field. This result can be also 
understood from a diagrammatic point of view in the following way. One first verifies that the 
expectation value of any observable can be expanded as a double series in G and N in which the 
power of N is always smaller or equal to that of G. One then verifies that the semi-classical value 
of this observable corresponds to the resummed series containing all terms governed only by the 
one-point function (T^). Moreover, all these graphs are 1-particle reducible and their weight 
is a positive power of GN . This analysis furnishes an alternative proof that the semi-classical 
treatment is the leading contribution and that it corresponds to a mean field treatment, see the 
Appendix for more detail. 

This diagrammatic analysis naturally leads to inquire about the next series. Upon having first 
summed up the leading series in powers of GN, one encounters a next series containing positive 
powers of G 2 N. This second series is governed by the two-point function (T^ u (x) T a p(x'))c, the 
'specific heat' of the radiation field. The second reason of having chosen the large N limit is 
that non-perturbative effects can be obtained by resumming this second series [1 9 . 

The third reason which has led us to choose a treatment based on a statistical basis arises 
from the trans-Planckian problem. Indeed, as previously discussed, the unbounded growth of 
frequencies encountered near the event horizon seems to invalidate perturbative treatments of 
radiative corrections, see [13] for an exemple of the divergences encountered. 

Let us briefly explain what is the nature of the trans-Planckian problem and why radia- 
tive corrections induced by quantum gravity might give some important effects when applied 
to Hawking radiation (or to quantum effects induced by the presence of an event horizon). 
When studying the origin of Hawking quanta one faces a difficulty which is specific to horizon 
physics: The configurations giving rise to Hawking quanta are characterized by ultra-high fre- 
quencies when measured by infalling observers near the horizon [301 EJ \12\ Hj. Indeed, in the 
semi-classical treatment, the use of free fields propagating in a classical background leads to 
unbounded frequencies as a direct consequence of the structure of the outgoing null geodesies 
near the horizon. Therefore, any ultra-violet scale which would signal the breakdown of the 
semi-classical treatment will be inevitably reached. 

This reasoning is nicely illustrated by considering sound propagation in an acoustic geometry 
which possesses a horizon: One first finds that the propagation is dramatically modified with 
respect to the standard propagation (governed by the d'Alembertian) when the FF frequency 
Q reaches £l c , the UV scale (the inverse inter-atomic distance). In particular, for all non- 
linear dispersion relations, one finds that the focussing effectively stops when this new scale is 
reached. Secondly, this dramatic modification of the near horizon propagation leaves no imprint 
on the asymptotic properties of Hawking radiation when the inverse inter-atomic distance is well 
separated [22j [23] from k, the surface gravity of the hole. 

Our aim is to show that similar results are obtained when computing non-perturbatively 
the gravitational effects driven by the connected two-point function (T^ T a p)c- We shall find 
that the trans-Planckian correlations which existed in the semi-classical treatment are washed 
out when the r — 2M reaches a, the new length scale which plays the role of the inter-atomic 
distance. Moreover, as already mentioned, this washing-out mechanism leaves the asymptotic 
properties of Hawking radiation unaffected: the thermal flux receives corrections which scale 
like (/-ccr) 2 and which are therefore negligible for large black holes. 
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1.2 The lesson from acoustic black hole physics. 

For the interested reader, we further discuss the relationships between our approach and the 
physics of sonic black holes. As just explained, the appealing feature of these models is to 
provide both a simple explanation (in terms of adiabaticity which essentially follows from scale 
separation £l c 3> k) for the robustness of the asymptotic properties of the flux, and a simple 
physical reason (a modified dispersion relation) which eradicates the ultra- high frequencies. (It 
should be pointed out that a similar trans-Planckian problem arises in inflationary models when 
studying the origin of the spectrum of primordial energy density fluctuations |25[ l32| [33] . In that 
case as well, scale separation and regularity of the metric are sufficient conditions to guarantee 
that the properties of the spectrum are unmodified [M].) 

Besides the robustness of the IR properties, the main outcome of these considerations is that 
a new universality has emerged: for all dispersion relations but the linear one, the blue shifting 
effect stops. Therefore, the never ending blue shifting effect obtained by using the linear (scale- 
less and non-dispersive) relation £1 = p now appears as an isolated and unstable behaviour. 
Thus, instead of asking: 

is Hawking radiation robust against modifying the dispersion relation ? 
we are led to question [35]: 

is O = p robust against radiative corrections ? 
or it is simply an artifact of free field theory ? 

These considerations suggest that quantum gravity should engender a new UV scale when 
evaluating radiative corrections in a black hole geometry. This new scale would then break 
the scale invariance of free field propagation and prevent the appearance of trans-Planckian 
frequencies. To verify this conjecture, one must determine the physical effects induced by the 
non-linearities engendered by gravitational interactions. When this is done, one can make con- 
tact with sound wave physics [37] by determining how phenomenologically describe by an effective 
linear equation (i.e., a non-trivial dispersion relation) the dissipative/dispersive effects induced 
by non-linearities. 

In this paper we shall implement the second question by computing how the radiative correc- 
tions driven by (T VV T VV ) modify the two-point function of out-going configurations. We shall see 
that O = p is not robust in that it provides a good approximation of the effective propagation 
law only when the FF frequency f2 obeys fi< 1/ct where a is given in eq. . 
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2 The Model 



For simplicity, we shall consider only s- waves propagating in spherically symmetric space times. 
For definiteness, the background metric is taken to be that resulting from the collapse of a null 
shell of mass Mq which propagates along the null ray v = 0. Inside the shell, for v < 0, the 
geometry is Minkowski and described by 

ds 2 = -(1 - )dv 2 + 2dvdr + r 2 (d9 2 + sin 2 ## 2 ) . (5) 

with M = 0. Outside the shell, the metric is Schwarzschild and given by eq. (|5]) with M(v) = M$. 
As we shall see, this choice of the collapsing metric will have no influence in what follows since 
the near horizon vacuum interactions are stationary. 



singularity 




Figure 1: Penrose diagram of the background geometry. The light-like infalling shell propagates along 
v = 0. The other continuous line emerging from J~ is v — Vh, the radial light ray which forms the event horizon 
after having bounced on r — 0. The dashed line represents a characteristic of the configurations <j>- which are 
responsible for Hawking quanta. The dotted line represents a partner's characteristic. For quanta reaching J + at 
late u, both of these characteristics are extremely close to vh, see eq. ©. The configurations (j>+ have support for 
v > and are always infalling. In this paper, we study the interactions between 0_ and cf>+ which occur outside 
the star, in the near horizon empty region, when the (initial) state of 4>+ and <f>- ls vacuum. 

To identify the degrees of freedom involved in these interactions, we first analyze the global 
properties of radial modes, when working in the geometric optic approximation, i.e. when 
working with d u d v <j) = 0. (In the exact d'Alembertian, see eq. ()10p . there is a potential around 
r = 3M which induces partial reflection, a phenomenon irrelevant for our purposes.) The ingoing 
massless waves fall into two classes according to their support on J~ , the light-like past infinity. 
The waves in the first class have support only for v < 0, inside the shell, and will be noted 
(f)-. They propagate inward in the flat geometry till r = where they bounce off and become 
outgoing configurations, see Figure 1. The relationship between the value of u of the geodesic 
which originates from v on J~ is [18j: 

V(u) = -4M(1 + e~ ku ) , (6) 
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The first class is thus divided in two subsectors: For v < —AM, the reflected waves cross the in- 
falling shell with r > 2M and reach the asymptotic region whereas those for > v > —AM cross 
the shell with r < 2M and propagate in the trapped region till the singularity. The separating 
light ray vh = —AM becomes the future horizon u = oo after bouncing off at r = 0. The 
configurations which form the second class live outside the shell, have support only for v > 
and are noted 0+. They propagate in the static Schwarzschild geometry, are always in-falling 
and cross the horizon towards the singularity. 

In Hawking's derivation of black hole radiation, the field operator obeys the d'Alembert 
equation. Hence the above classical properties apply: The configurations for v < v h give rise 
to the asymptotic quanta, those for v h < v < to their partners [5] whereas the configurations 
described by 0+ play no role in the asymptotic radiation. This follows from the asymptotic 
(ku 3> 1) behaviour of the relation V(u) 

V{u) -v H oc e~ KU . (7) 

As shown in [1], this exponential is responsible for the thermal radiation at temperature k/2it. 
It also shows that Hawking quanta emerge from trans-Planckian frequencies on J~ since codV = 
Xdu (where u> = id v on J~ and A = id u on J + ) gives u> oc Xe KU . Finally it fixes the correlations 
between the asymptotic quanta and their partners, see eq. ([2]). These follow from the fact that, 
on J~ and in the vacuum, the rescaled field (f> = V Airr 2 x (where x 1S the 4D s-wave) satisfies 

/•OO A, . 1 

{(j>(v) cf>(v')) = / ^L e -Mv-V) = _±_ ln r v _ y f_ ie) + constant . (g) 
Jo Airto An 

Since this equation is valid for all v,v' one might think that there also exist correlations 
between </>_ and (/>+. However, for late Hawking quanta, they effectively vanish since these 
quanta and their partners emerge from configurations which are localized extremely close to vjj 
as indicated in eq. ([7]). For a description of the other properties of Hawking radiation, we refer 
to the review [3]. 

In brief, in the absence of gravitational interactions, </>_ and (/>+ are effectively two inde- 
pendent fields. By independent we mean that by sending quanta described by wave packets 
built with there is no induced emission of Hawking quanta. Indeed, in order to get induced 
emission [31] at time u, one should send </>_ quanta localized close to vjj as indicated in eq. ([7]) 
and correspondingly characterized by high frequencies u oc Xe KU ^> k. 

Let us now analyze more closely how these properties translate in Fock space. When evalu- 
ated in the background eq. ([5]), the action of 4> is 



Sf = - Idvdr 



(9) 



with M(v) = for v < and M{v ) = Mq for v > 0. Being interested in the near horizon physics, 
we have dropped the potential term of s-waves, (2Mo/r 3 )0 2 , since it does not affect the near 
horizon propagation. This can be seen by using the double null coordinate system u = v — 2r*,v. 
Using them, the 4D d'Alembertian reads 



(10) 



where 4>i is the rescaled mode of angular momentum /. Thus, as emphasized in p3], the propa- 
gation of waves (at fixed angular momentum and even for an arbitrary mass) effectively obeys a 
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2D conformal invariance in the near horizon geometryLO This is confirmed by the fact that, clas- 
sically, the trace of 2D part of vanishes independently of the equations of motion. Thus, in 
our model for s-waves, has only two g-number components, T vv = (d v (j)) 2 and T uu = (d u (j)) 2 . 

The 2D conformal invariance also implies that the Fock space is composed of tensorial prod- 
ucts of in-falling states (on which (f) + acts) and outgoing states. In a Schroedinger language this 
means that an initially factorized state (i.e., = |^+) ® |^-)) remains factorized when its 
evolution is governed by eq. ([9]). This factorizability explains the absence of induced emission 
when adding (f> + quanta. In a Heisenberg language, it tells us that any matrix element of <j) is a 
combination of matrix elements of </>_ and (f>+ evaluated separately. This implies in particular 
that the connected part of the two-point correlation eq. (pQ) identically vanishes for all factorized 
states. This applies to the "Unruh" vacuum, the state describing Hawking radiation. Physically, 
the vanishing of eq. ([I]) means that the fluctuations of T vv and T uu around their mean values 
are completely uncorrelated. It should be clear that this absence of quantum correlations is 
precisely what is contested by t'Hooft, see footnote 3. 

Finally, we mention that, in spite of this absence of correlations, the mean value of T vv 
and T uu are related to each other by energy conservation through the 2D trace anomaly |36j . 
However, this third component of T^ v does not fluctuate: it is a c-number. Hence it cannot play 
any role in the gravitational interactions between the field operators </>_ and <\>+. 



3 The gravitational interactions between and cj) + 

The aim of this Section is to obtain the dominant part of the action governing the gravitational 
interactions between </>_ and </>+. In the next Sections, we shall study the consequences of these 
interactions with particular emphasis on the correlations they induce. 
The generating functional governing our matter-gravity system is 

Z = JvcpVh e l[S t+h +Sh '° ] . (11) 

h is the change of the metric with respect to the background g discussed above and Sh, g is the 
Einstein-Hilbert action. (We only consider metric fluctuations in the Ricci flat region outside 
the infalling matter forming the hole.) S^ +h is the action of <f> propagating in the fluctuating 
geometry g + h. 

When the metric fluctuations are spherically symmetric, h can be characterized by two 
functions: ip, fi. Moreover both are completely determined by the energy-momentum tensor of 
cf). This is like the longitudinal part of the electro-magnetic field which is constrained to follow 
charge density fluctuations, by Gauss' law. The line element in the fluctuating metric can be 

7 This invariance leads directly to the trans-Planckian problem: the steady production rate of outgoing quanta 
arises from an integral over in-frequencies to whose measure is that of a 2D massless field. Explicitly one obtains 
that the thermal distribution is multiplied by dui/u — ndu since u tx e KU , see [3l 113]. 
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written as [18J 1 



i 2 V 



1 



2M 



dv 2 + 2du dr 



+ r 2 d^ . 



dv dr 



d v (f>d r (f) + 1 



2M\ {d r <f>f 



(13) 



where M = Mq + r) for v > 0. 

In this new metric, the matter action is the same as in equation eq. (|9|): 



(14) 



The new mass function M incorporates the only relevant metric change /x. Indeed, Sg +h is 
independent of ip, thereby recovering the 2D conformal invariance mentioned earlier. 

Our aim is to work out the first order corrections due to the gravitational interactions 



between 



and 



To this end only quadratic terms in h should be kept in Sh «• The 



Gaussian integration over h can be performed (this is equivalent to solve the linearized Einstein's 
equations). It gives rise to a self-interacting field theory described by 



Z 



V(j)e 



isf -\-iS{nt 



(15) 



By construction Sfat is a non-local quadratic formal of the energy-momentum tensor of <fi. 

To identify the relevant part of S{ n t, we first recall that T^ u has only two fluctuating compo- 
nents, thanks to the 2D conformal invariance. Thus, in a perturbative treatment (such as the 
interacting picture) one obtains two types of interaction terms only. First one has self-interaction 
terms depending on 0_ or 4> + separately. These terms do not destroy the factorisability of the 
theory and will not be considered in what follows F^l Secondly, one has a cross term coupling 
0_ to (f> + . This term will inevitably break the factorisability of the theory into the ± sectors. 
Therefore, the two-point function (pQ) will no longer vanish. Let us analyze this term in more 
detail. 



This line element differs from that used by Bardeen [6]: 



ds 2 = 



l- 2Mo+2 ^W + 2d,dr 



+ r 2 dQ.l 



The ip function is the same whereas, to first order in ip and fiB, the mass fluctuation n — hb — ip(r — 2Mo)/2. 
The usefulness of our choice is that ip no longer affects the null geodesies. When expressing T M „ in terms of the 
null fluxes T VV ,T UU , Einstein's equations give 



d r **p = 4T uu /(r-2M). 



(12) 



9 In the t,r coordinate system, i.e., when g r t = 0, Sint is given by a linearized version (see equation (90) in [5]) 
of the BCMN [38] Hamiltonian. 

10 The validity of this simplification (also adopted in [121 1141 [T5] l requires further analysis. On-shell, the tf>+ <f>+ 
contribution to Sint vanishes. This can be understood from the fact that the Vaidya metric (2) is an exact solution 
for any classical infalling massless flux T vv (v). The <f>- (f>- contribution to Sint is more tricky to handle in the 
advanced coordinates v,r. The reason is that infalling geodesies are affected by the presence of an outgoing flux 
T uu (as clearly seen when using the coordinates u, r) . This modification translates in v, r into a deformation of the 
description of outgoing geodesies u = u(v,r) and it is this effect which is responsible for the ^>_^>_ contribution 
to Sint- Let us finally notice that a non-perturbative treatment of the self-interactions of <f>- has been developed 
in [39]. It leads to small effects 0(k/M) and induces no damping of the waves when approaching the horizon. 
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Since infalling configurations obey <9 r </>+ = even in the presence of gravitational interac- 
tions, the cross term coupling </>_ to <p + is, see eq. (UU), 



Smt = G [°°dr [°°dv ^(d r ^f 
Jo Jo r 



oo roc 



JO 



(i+(y) du 2 



G dr dv^^(—YT uu . (16) 



r dr 



where G is Newton's constant. We have re-introduced it in the front of ^ to read more easily the 
order of the interactions between 0_ and <j) + in the forthcoming equations. /j.+ (v) is the mass 
fluctuation driven by 4>+- Einstein's equations constrain is to be 

H+(v) = f^v' T vv (v') = f dv' (d v ,<t> + ) 2 . (17) 

JO JO 

(The reader might be surprised by the fact that we are using on-shell fields <f>± in Sint, i.e. that 
we have used the equations of motion. In principle indeed, only the off-shell field cf) should be 
used in the action. However, when calculating perturbatively lowest order corrections in G, this 
amounts to use eq. (fT6j) as it stands.) 

We are now in position to show that the gravitational interactions between ([>+ and 4>- 
diverge on the horizon. To this end, let us consider two classical fluxes described respectively by 
T vv = uj 5{v — vq) and T uu = A S(u — no), ui and A are the asymptotic energies measured on J~ 
and J + respectively, and vq and uq are such that the two spherical shells meet at ro in the near 
horizon geometry, for ro — 2M <C 2M . In this case, using ro — 2M ~ 2A4e K ^ v °~ u °\ one obtains 

S int - 4G — - . (18) 
r /2M-l 1 1 

The action S{ n t diverges as ro — * 2M like the FF frequency f2 did in (3). The difference with 
(3) is that Si n t is a scalar. Hence the divergence in eq. (fT8|) is coordinate (gauge) invariant. 

To get an estimate of where the gravitational interactions become strong, i.e., can no longer 
be ignored, let us consider two shells whose asymptotic energy is Hawking temperature, i.e. 
to = A = k. The condition Sint = 1 is reached for 

r /2M - 1 = AGk 2 . (19) 

The proper distance (ro — 2M ~ 1/M) is much smaller than Planck length, see footnote [U 
This simple estimate will be recovered in Section 6 when considering radiative corrections in the 
vacuum. We should emphasize this last point: even though our approach closely follows that 
of [12] (it can be considered as an s-wave reduction of it) we shall not study the interactions 
between (j> + and (ft- quanta. Rather we shall focus on the residual interactions when the state 
of </>+ is vacuum. For earlier attempts in this direction, we refer to [141 115j . Before analyzing 
these second quantized effects, it is instructive to solve two preparatory exercises with on-shell 
fluxes. 

In the first we shall show that Sint ac ts only as a shift operator of the asymptotic value of 
u. In spite of this simplicity, in the second exercise, we show that Si n t nevertheless engenders 
an entanglement which prevents the factorisability of the states into db sectors. This provides 
an explicite exemple of a quantum effect induced by eq. (|16p which cannot be described in the 
semi-classical treatment. 
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4 Non-vacuum gravitational effects 



4.1 Classical interactions and shifts in u 

We first consider the following problem: Given two classical field configurations: $L specified 
on J~ for v > and cjp_ on J + , what is the value of the field amplitude (j) near the horizon 
when taking into account the gravitational interactions of eq. (|16p ? 

Because of the 2D conformal invariance, (j) still decomposes as 0+ + </>_. Then, since 9 r 0+ = 
is exact in our gauge wherein v stays light-like, 4>+(v) = <fi+(v) to all orders in G. Thus the 
infalling flux of energy is unaffected by the energy carried by <f>— and it is given by its initial 
value on J~: T vv = (<9„</ , +) 2 • Hence, u + of eq. (fTT|) acts as a given metric change in the equation 
of motion of </>_: 

2tf l + 2c^; 



_ = . (20) 
Since this equation is linear in 0_ and first order in the space-time derivatives, its exact solution 



is 



Mv,r) = 0° r)), ( 21 ) 
where u^(u,r) gives (in the coordinate system of eq. (fT3"])) the outgoing null geodesic in the 
modified metric characterized by Mq + Gfi + (v). 

To determine the function u^{v,r), one notices that it also obeys eq. (|20|) with the (final) 
boundary condition that it converges to the un-modified value uo(v,r) = v — 2r for r — * oo. 
Hence, to first order in G, the change <5u = u^—uo is determined by a non- homogeneous equation. 
Using the fact that 2d v + (1 — 2Mo/r)<9 r defines 2d v \ Uo (by definition of uq(v, r) = constant), 5u 
obeyspEE] 

d v \ U0 5u = G — d r \ v u . (22) 
r 

The solution is 

5u(v)\ uo = G Hdv' 1^ +{V '\ M , (23) 
Jv r(v')\ Uo -2M 

where r{v)\ Uo is obtained by inverting uq(v,t) = v — 2r* . As in the action (|16|) the integral in 
eq. (|23p is dominated by the near horizon region when the denominator r — 2Mq <C 2Mq. A 
tiny infalling energy flux /x+ <C M can therefore induce an arbitrary large change of u. 

The lesson we got from eq. (|2ip is that the eikonal approximation is exact: the scattered value 
of the field amplitude is given by its asymptotic value evaluated along the modified characteristic 
u^(v,r). Thus, classically, the gravitational interactions encoded in eq. (|16|) only induce a shift 
of the argument of field. They do not induce non trivial non-linearities in the field amplitude in 
that fi+(v) can be computed irrespectively of the value of 0_. The origin of this miracle is the 
2D conformal invariance. 

4.2 Entanglement between 0_ and <p + 

In spite of this absence of non-linearities in field amplitude, we shall now prove that the quantum 
evolution governed by the action + S% n t dynamically engenders entanglement between the 
otherwise uncorrelated ± sectors. To this end, we consider the evolution of an initially factorized 
wave function 

= |*») <g> |*~) . (24) 
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The infalling part \^f l J^) is specified on J for v > 0. To clearly exhibit the entanglement, we 
choose to be a superposition of two well defined and separated states: 

\V™) = A\y™' a ) + B\V™' b ). (25) 

The two kets are normalized and orthognal to each other: (fy 1 ™' 1 ^ 1 ™' 3 ) = 5^. Thus A,B are 
probability amplitudes obeying \A\ 2 + |i3| 2 = 1. By well defined and separated we mean that 
the two fluxes associated with each component, (T£ v ) = ' l |T OT |\l/!, ' ) with i = a,b, are well 
localized in v and separated from each other. 

The other piece of the initial ket, l^!! 1 ), is specified on v < 0. After reflection on r = 0, it 
determines the state of outgoing configuarations. For the moment we do not need to further 
specify its quantum state. 

Having specified the initial state, we study the quantum dynamics (in the interacting picture). 
The evolution operator e lSint acting on the initial state |v]/ m ) gives two uncorrelated evolutions 
weighted by A and B: 



|\J/\ = giSint |xj/«n\ _ e iSi„t 



(A\V™' a ) + B\¥? b ))®\i!> i 



Ae l 



+ Be l 



(26) 



One should thus study each piece separately. This is nothing but the expression of the superpo- 
sition principle. 

In each state, the equation <9 r 0+ = 0, now viewed as an Heisenberg equation, tells us that 
the evolution in the + sector is trivial, as in classical terms. To study the evolution of the 
— sector, we perform the approximation which consists in neglecting the fluctuations of T vv 
in each infalling state I^V 1 ''). In this approximation, the evolution operator e lSint becomes a 
c-number for infalling configurations and acts only on outgoing configurations. One thus has 



A (V + "< a > ® e^H*™)) + B (\*™' b ) ® e iS int\V™)) . (27) 



In the ^4-weighted state, the operator e lS ^ governs the propagation of 4>- in the o-modified 
metric characterized by the mass function M a = Mq + Gfj,i_ with 



V 

in, a 



M»=/ dv'(^ a \T vv (v')\^ a ), (28) 
Jo 

whereas in the B state, one finds the 6-modified metric characterized by fA_(v). The evolution 
in each case is thus governed by the time-ordered product of "its" evolution operator acting on 
same initial l^!! 1 ). Explicitly, the two hamiltonians acting on |^™) are given by eq. (|16p with 
the corresponding the c-number metric changes fj,\_ , i = a,b. In the above treatment, we took 
into account the first moment of T vv , {T vv } 1 , the mean value in each state. We indeed neglected 
the highers moments which govern the fluctuations in each state. In doing so, we accounted for 
the fluctuations of T vv which are due to the fact that the infalling state is a superposition, but 
only those. This will be clarified below. 

The entanglement induced by Si n t acts, as usual, as a measurement. For the interested reader 
we recommend the very instructive reading of Chapter IV. in [40 1. Consider, as in that chapter, 
the Stern-Gerlach experiment wherein the center of mass motion of an electron in a magnetic 
field is governed by its spin projection along that field. The mapping from that situation to 
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the present one is as follows. The two kets representing the spin projections are here played 
by the two infalling states ). The center of mass wave function is played by the outgoing 
wave function ) and the interaction hamiltonian is Si n t of eq. (|16p . The analogy works 
quite well when the initial outgoing wave function |<J/™) also describes a localized flux. Then, 
its "image" on J + would be either a spot at uq + 5u a with probability \A\ 2 , or one at uq + Sub 
with probability 1 — \A\ 2 . The location uq is that of the spot when the gravitational interactions 
are ignored, whereas the values of the shifts 5u a ,Sub are given by eq. ([23]) . evaluated along the 
unperturbed geodesic uq, fed with the mass changes /x+ and /i b + respectively, and with the lower 
value v set to 0. 

This quantum result should be compared with what would have been obtained by applying 
the semi-classical treatment to the entire wave function (rather than to each sector separately). 
In that case, the change in the common location is driven the mean mass change 

J, + {v) = \A\ 2 ^ + {v) + \B\ 2 ^ b + {v). (29) 

Therefore, the semi-classical treatment incorrectly predicts a single spot on J + located at the 
"mean" position uq + 5u with 5u = \A\ 2 5u a + \B\ 2 5ub- 

In the above exemple, the validity of the semi-classical treatment rests on the possibility of 
neglecting the fluctuations of the <^+-operator-valued shift 5u. The dispersion about the mean 
quantifies their importance. One finds 

' «o) 

2 



(5u UQ 5u U0 ) c = ((6u Uo - {5u Uo )) ) = {(5u U0 ) 2 ) - (5u % 

= \AB\ 2 (5u a (u ) - 5u a (u Q )f , (30) 



= AG 2 rr^dvt ,, „ ^t'f/ft TFT ■ (31) 

Jo Jo (r(tfi)k, -2M ){r(v 2 )\ U0 - 2M ) V ' 

To get the second line we have used the results of the above two paragraphs. From eq. (|30p one 
verifies that the semi-classical treatment is valid iff 

{^Ju„)a_ IABl A e "-- s ";)\^ (32) 



OC 



It is negligeable either when |^4i?| 2 — > 0, that is, when one of the infalling flux is rarely found, 
or when the two shifts in u are close to their mean 1^1 

To get eq. ([3T]) we simply have used the definition of 5u, eq. ([23]) . and that of connected 
correlation functions, see ([1]). Eq. (131|) displays the relation bewteen the dispersion of 5u and 
(fi + (v)fj 1+ (v')) c , the correlation function of the metric fluctuation fi(v) in the quantum state 
one is dealing with, here given in eq. (|25p F^l Because of the horizon, the importance of the 



11 We refer to 56, 57, 58 for a treatment of the solutions of the Wheleer-deWitt equation in which the gravita- 
tional back-reaction is computed for each component of the matter wave function separately. This radically differs 
from the usual mean field treatment in which the gravitational back-reaction is determined at once[59]. Only the 
first treatment respects the supersposition principle. Indeed the back-reaction in the second case contains the 
probabilities to find each component, exactly like in Su = |y4| 2 5u a + \B\ 2 Sut,. Hence the limitations of the validity 
of the second (the semi-classical treatment) are not intrinsic to the gravitational dynamics. 

12 We remind the reader that to get (|30[) we have neglected the fluctuations of fi in each infalling state, i.e. only 
the dispersion induced by the A-B superposition has so far been taken into account. The aim of the next Section 
is to incorporate the inherent fluctuations which are present in any infalling state including the vacuum. This 
explains the title of the present Section: Non-vacuum gravitational effects. 
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dispersion is determined by two things. On the one hand, the fluctuations of 5u are driven by 
the fluctuations of the mass function p, = J dvT vv which are themselves determined by those 
of the infalling flux T m , given the state of infalling configurations. The connected correlation 
function (T VV T VV ) C quantifies their dispersion. On the other hand, the mass fluctuations are am- 
plified by the denominators in eq. (f3Tj) . Therefore, even if the relative importance of (T VV T VV ) C 
with respect to the square mean (T vv ) 2 is small, the semi-classical treatment can give incorrect 
predictions when the fluctuations are sufficiently amplified by the denominators. Instead, in the 
absence of amplification, the mean theory provides a reliable approximation, unless of course if 
Ma(^) ~ /J-biv^/M ~ 1, but in this case the a,b branches completely decouple/decohere, and one 
should consider each possibility separately. 

In brief, the crucial points we have reached are the following. First, unlike what one encoun- 
ters in usual cicumstances, i.e., without an event horizon, tiny fluctuations of T vv might give 
rise to large shifts in u because of the amplification due to the growth of the gravitational inter- 
actions when the configurations meet near an event horizon. Second, these fluctuations break 
the factorizability of the theory into the (f>+,(f>- sectors. Third, the quantity which governs the 
validity of the semi-classical scenario is {T VV T VV ) C , the connected two point function of T vv . 

The aim of Section 5 is to extend the analysis of eq. (|3ip when the infalling configurations are 
in the vacuum, i.e. in a coherent state as opposed to a superposition as in (|25p . In that Section, 
the (inherent, minimal, and divergent) vacuum fluctuations will be the only contribution to 

4.3 Relationship with former treatments 

Before considering vacuum effects, it is also interesting to relate eq. f)15|) to the former treatments 
of black hole evaporation discussed in the literature: Hawking's approach pQ and the semi- 
classical treatment. 

Hawking's approach formulated in a fixed geometry is recovered by putting G[i+ = in 
eq. ()15p . Then Z factorizes as Z + ® Z_ (when ignoring the trace anomaly) and 0_ is a free 
outgoing field propagating in the background geometry g. Thus </>+ drops out from all matrix 
elements built with the operator </>_. It should be emphasized that the trans-Planckian problem 
encountered in Hawking's approach directly follows from this factorisability. Indeed it is the 
absence of gravitational coupling between the + and — sectors which permits the unbounded 
growth of frequencies when probing, near the horizon, configurations specified on .T + [5| [18]. 

The semi-classical treatment [HI El IE] can be obtained from the path integral formalism eq. 
(jlip by first integrating over at fixed h and then searching for the classical extremum of h. 
In this approach, by construction, the fluctuations of h and are neglected. Thus the near 
horizon propagation of 4> is governed by a single but now self- consistent metric governed by 
mean (//_)_ (v)}. This mean evolution characterizes by the shrinking of the horizon area according 
to 



dv 



(T V v)\r=r horizon =2M • (33) 



When working in the vacuum eq. ©, the (properly subtracted [3]) expectation of T m is 

(T m (v))\ r=2M{v) =-^(^^j , (34) 

where k(v) = l/4M(y) with M(v) = Mq + G{fi+(v)). This flux has the opposite value of a 
2D thermal flux. The only change with respect to the fixed background approach of Hawking 
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is the replacement of Mq by M(v). Thus the propagation of out-going configurations is hardly 
affected by the evaporation as long as it is slow, i.e., as long as M(v) 3> Mpi anc f.. 

Therefore, in the semi-classical scenario, the trans-Planckian problem stays as in Hawking's 
approach: The coupling between 0_ and the mean change (fJ>+) is incapable to provide a taming 
mechanism since it does not open new interacting channels. To solve this problem clearly requires 
to take into account the fluctuating character of the interactions between (/>_ and i.e., the 
possibility of entangling their wave functions, as in the quantum mechanical exercise presented 
above. 

5 Modified two-point function 

Our aim is to determine how the (dressed) two-point function of <f>- 



where Z is given in eq. (|15p . is affected by the gravitational interactions encoded in Si n t when 
the infalling configurations are in their vacuum state. 

To evaluate eq. (j35[) beyond the semi-classical treatment, one should adopt some rules 
to cope with the UV divergences. Even at one loop, there are several inequivalent graphs 
which result from the various Wick contractions. To extract consistently the contribution of the 
simplest ones, we propose to consider N copies of </>. The calculation of eq. (I35h can then be 
achieved in two different (but equivalent) approaches. The first consists integrating first over 
the TV — 1 ~ N spectator copies not appearing in the numerator in so as to determine the 
(linear term in N of the) influence functional (IF) [UJ governing the effective dynamics of <fi. 
The other consists in developing e int in both integrands of eq. (135j) in powers of Si n t so as to 
engender the (connected) graphs governing the radiative corrections. Then in the large- N limit 
the graphs can be classified into infinite series according to the relative power between G and 
N, see the Appendix for the details. The leading series (in GN) reproduces the semi-classical 
treatement. The next series of graphs, those weighted by powers of G 2 N, all correspond to the 
graphs engendered by the term in the influence functional which is linear in G 2 N. 

The usefulness of the IF approach is that non- linear (in G 2 N) modifications of the propaga- 
tion of outgoing configurations are easily taken into account through this influence functional. 
The same results can of course be reached from the diagrams approach at the cost of re-summing 
the corresponding infinite subsets of graphs. It is in the identification of these infinite subsets 
that the large N limit finds its justification. For a schematic description of these diagrammatic 
aspects, we refer to the Appendix. In what follows, we shall pursue with the IF approach. 

When computing the lowest order corrections to the effective action in eq. (135H . we can use 
eq. (|8l), the 'free' propagator of </>+El This approximation concerning degrees of freedom not 
directly involved in the matrix elements (i.e., which factorized out in the absence of interactions) 
is a common procedure both in quantum field theory where it gives the vacuum contribution, 
see chapter 9 in |41j . and in statistical mechanics (e.g., the polaron, chapter 11). In our case, in 
this approximation, the IF gives rise to a non-local action which is a sum of terms containing 

13 In would also be interesting to determine if the higher order terms in powers of G will become operator- 
valued [14| in (/>_, thereby obtaining a situation analogous to that of transition amplitudes when enlarging the 
phase space so as to take into account recoil effects [47] . 




(35) 
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(d r (j)-) 2 and kernels given by the Wick contractions of T vv evaluated with eq. ([8]). To order 
G 2 N, one obtains^ 

Sif = i(S in tS in t) + = iG 2 N fd 2 x [d 2 x'(rr')-\d r <f>-) 2 (»+{v)n+{i/)) (cV0_) 2 , (36) 



where ( )+ means that the expectation value applies to 0+ only. Using eq. ([8]), the connected 
two-point function is 

(T vv (v) T vv {v')) c = J— . (37) 



Then, eq. (fT7|) gives 



(fi+(v) fj,+ (v')) 



1 1 



96-7T 2 (v — v' — ze) 2 
1 Z" 00 

^2 / dm uj exp[-iuj(v - v)] . (38) 



o 



This equation gives the mean fluctuations driven by one 0+ field in the unperturbed {G = 0) 
vacuum state, see |17| for a analysis of the 4D two-point function of induced metric fluctuations 
in Minkowski vacuum, see also [H] for a analysis of metric fluctuations in the near horizon 
region. 

It should be noticed that ( / u + (w)/i + (w / )) is not real. This follows from the quantum vacuum 
which contains only positive frequencies when hit by <f>, see eq. ©. Notice that eq. (|37p gives the 
"un-subtracted" value of the connected two-point function. As shown in [43], the counter-terms 
which lead to the renormalized one-point function, eq. (|34p . provide divergent contributions to 
eq. (|37p which tame its singular behaviour as v — > v' , and make it a well-defined distribution [17]. 

Keeping only eq. (|36p in the IF (or by summing the corresponding infinite set of Feynman 
graphs, see the Appendix) is equivalent to work with a stochastic (i.e., a classically given, 
albeit not real) and Gaussian ensemble of metric fluctuations. By equivalent we mean that 
the graphs obtained from the Gaussian integration in the stochastic treatment are in one to 
one correspondance of those of the £r 2 iV-truncated quantum treatment. Therefore all matrix 
elements of such as eq. (|35p . can be computed from the stochastic theory. In brief, as far as 
the propagation of 0_ is concerned, neglecting the time ordering in the evolution operator, and 
to lowest order in G, the functional integration over <p + in eq. (|35p defines an effective stochastic 
ensemble of metric fluctuations governed by eq. (i38|) . 

In this case, all the techniques developed in [18J apply. In what follows we shall present 
schematically the main results and we refer to this work for details. The key point is the 
following. Because of the Gaussianity of the ensemble, one can obtain non-linear corrections to 
eq. (|35l) from the fluctuating characteristics of eq. (f20|) . i.e., the outgoing null geodesies u^(v, r), 
the non-trivial solutions of ds 2 = in the fluctuating metric eq. (|13p . In this we recovered that 



14 In the above expression, we did not take into account the time ordered character of the evolution operator 
e lSint . It would be interesting to determine the consequences of this neglect. We refer to App. A of I45l[ for a 
discussion of this point in the context of atomic transitions. We also refer to the recent work of Weinberg l4&j 
wherein the 1-loop radiative corrections to the two-point function determining the power spectrum in a inflationary 
context are computed (in a large N limit as well) taking properly into account the time ordering. In \21f . we tried 
to improve the calculation of S19j by taking it into account. We obtained an extra divergence in eq. which 
induces dispersive effects. However the meaning of this extra term is unclear as it seems to correspond to a tadpole 
contribution which should be subtracted. 
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there is no non-linearities in the field amplitude: as in Section 4 the non-linearities in G only 
occur through the characteristics. 

To determine the effects engendered these metric fluctuations, it is instructive to analyze 
the backward in time propagation of configurations representing asymptotic Hawking quanta. 
In particular it will reveal the role played by shift in u studied in section 4.1. To this end, we 
consider the Fourier transform, performed on S , of the in — in Green function (that obtained 
by taking the expectation value in the initial vacuum on J~ ): 

G-(X;v,r)= [ dm e iXui G-(y x = oo, ui; v, r) . (39) 



In the unperturbed metric, u = 0, i.e. ignoring 5j nt in eq. (|35j) . we get 

G-(X;v,r) oc e iXu(v ' r \ (40) 

where u(v, r) = Uq(v, r) = v — 2r* is the unperturbed null characteristic. Hence near the horizon, 
the Fourier transform of G_ behaves as 

e -i\u {v,r) _ e ~tXv ( r _ 2M ) iKX . (41) 

It possesses an infinite number of oscillations as r — > 2Mq with a radial momentum given by 

p r e ~^Mv,r) = _ idr e -»A«o(«,r) = e -iXu (v,r) _ / ^ e -i\u (v,r) _ /^x 

r /2Mq — 1 

A FF observer will attribute to this wave a FF frequency that grows as in eq. (|3|) because 
dr\ v oc dr, where r is his proper time. In fact p r should be interpreted as a frequency rather 
than a momentum (its sign fixes that of the Klein-Gordon inner product [23, 26J). 

Taking into account the first order change in u, see eq. (|23p . in any Gaussian ensemble of 
infalling metric fluctuations, the Fourier transform of the ensemble average of the in — in Green 
function of eq. (|39h will be governed by the averaged waves given by 

.2 

ll e -i>^iM{v,r) \\ _ e -i\u (v,r) & - ^- {{Su(v)Su(v))) _ 

In the above the averaging acts on the fluctuating quantity /u + (i>). 

We now focus on the particular (Gaussian) ensemble wherein the fluctuations are those 
induced by ./V fields in the vacuum, i.e. they are governed by iV times the expression of eq. ()38p . 
Then, using the fact that r (v)\ uo - 2M ~ 2M e K{v - u °\ one obtains 

G 2 N f A d W k 2 oj 

(du(v)\ uo du{v)\ uo ) 



(r/2M - l) 2 Jo 3 k 2 +lj 2 

= ^ (r/2M - If ■ (44) 
The spread a\ governs the damping of the backward propagated waves. It is equal to 



GktJn 1ii(A/k)/3, (45) 



when the hard UV cut-off A satisfies A>k. We have introduced A to define the two integrals 
over ijj. Notice that A is a Lorentz scalar in the sense that it is the energy of an s-wave in its 
rest frame in a stationary and spherically symmetric background. 
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The important result of eq. (|44p is that a\ is not proportional to A even though — A 2 . 
Notice indeed that a\ hardly depends on the value of A since 

OA=4M = V2 °A=Mpi anck ■ (46) 

This insensitivity follows from the fact that high frequencies (to ^> k) are damped by the 
integration over v' in eq. (|23p . The frequencies uj ~ n dominate the contribution to a in eq. 
(|44p . see [46J for a similar result in an inflationary context. However they are not sufficiently 
damped to give a finite result. Therefore, in the simplified treatment we are using, the value of 
A must be chosen from the outset0 Instead, in an improved treatment of the regularization 
of the divergences in eq. (144D . we believe that this ambiguity will be resolved. That is, when 
using the properly subtracted [33] two-point function {T^ u T a p), only uj ~ k should contribute 
to the (j's. The reason is the following: in the UV domain, all expressions become Minkowskian 
in character and hence cannot contribute to Lorentz breaking effects such as those engendered 
by a. 

The second result of eq. (|44p is that (5u Su) diverges as r — ► 2Mq. Thus the correlations 
between asymptotic quanta and early configurations, which existed in a given background as 
shown in eq. (I4ip . are washed out by the metric fluctuations once r — 2Mq ~ a\. The reason 
of this loss of coherence is that the state of 0_ becomes correlated to that of \12\ I14j . 
Physically, this loss of coherence implies that induced emission |31| n o longer exists when the 
threshold energy (measured in the FF frame) 1/a = 4Mn is reached] 16 ! Phenomenologically this 
loss can be viewed as a dissipation of outgoing waves, and, as in condensed matter [37] [16], it 
can be described by an effective dispersion relation. Explicitely, using the notations of [26[ [27]. 
we get an effective relation Q 2 = F 2 (p), where Q and p are respectively the frequency and the 
radial momentum measured in the FF frame, and where F 2 behaves as 

F 2 (p) =p 2 (l- 2i a 2 up) =p 2 (l- 2i — ) . (47) 

Pc 

To obtain this, it suffices to rewrite the r.h.s. of (|43p as e~ lXu ° e~°" 2p o/ 2 and to work in the 
momentum representation. Then with the help of eqs. (123, 124) in [26] the identification 
of the function F is straightforward. It should be noticed that the UV momentum scale p c 
which weighs the cubic term is p c = (cr 2 ^)^ 1 = (1/a) (Lpin)~ 2 = Mpi(LpiK,)~ 3 , and not 
p c = 1/a as one might have thought. This mismatch of UV scales illustrates that it is probably 
meaningless to try to identify a well defined dispersion relation when starting from the result: 
the modified propagation of eq. (03]) which unequivocally states that backward propagated 
modes are dissipated when their FF momentum reaches 1/a. The reason why this identification 
fails is that the modified propagation arises from the non-trivial properties of the background 
and not from short-distance physics as it is the case in condensed matter. 

We should further explain the physical relevance of these results. To this end, one must 
identify the matrix elements of <f>- which are sensitive to the metric fluctuations (and governed 
by the ensemble averaged waves eq. ()43p ) and those which are not. The simplest exemple of an 

15 In spite of this ambiguity, we still believe that the uj ~ k dominance and the associate robustness of a give 
credit to the validity of |^5p. This is conforted by the fact that a also governs the modifications of the asymptotic 
properties of Hawking radiation' 18 1.) 

16 An interesting and unsolved question raised by this loss is whether new correlations are induced by the 
gravitational interactions as the same time as the old ones are washed out. This phenomenon of replacement of 
correlations has been clearly derived in a slightly different context in 14 9f - 
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operator which is sensitive is the Fourier transform of the in-in two-point function in eq. (|39p . 
On one hand, the phase of the wave function used in the Fourier transform is not affected by 
the metric fluctuations. On the other hand, the location v, r of the second field operator is also 
unaffected by the metric fluctuations. However the two-point correlation function is sensitive to 
the metric fluctuations encountered from J + to v , r0 It is this (unusual, see below) discrepancy 
between the phase at each point which explains why the ensemble averaged waves of eq. fj43f) 
govern the above Fourier transform of the two-point function. 

Instead usual expectation values, such as for instance the in-in Green function with two 
points evaluated at fixed u on J + (or two nearby points close to the horizon), are not severely 
affected by the metric fluctuations because these expectation values are not governed by the 
ensemble averaged waves (I43D . The reason is that the ensemble average is performed after 
having computed the (ft- expectation value for each member of the ensemble. (This is not a 
choice: the stochastic ensemble is merely a tool to reproduce quantum mechanical expectation 
values. Its quantum origin fixes the rules of ensemble averaging without ambiguity). In our case, 
this implies that the shift eq. (123j) affects coherently the phase at each point, see Section IV. A in 
|18j . This guarantees that the shift drops out in the coincidence point limit. This cancellation in 
turn guarantees that the asymptotic properties are (almost) unaffected since the Green function 
possesses the usual Hadamard singularity [12]. By almost we mean that the corrections scale 
like (kct) 2 and thus are order 1/M 4 . It is important to point out that it is again the dynamically 
induced scale a\ and not the UV cut-off A which governs these corrections. 

We would like to further discuss the fact that the metric fluctuations strongly affect the 
correlations between configurations specified on J + and near the horizon without modifying 
the short distance behaviour of the Green functional The radical difference of the impact of 
vacuum gravitational interactions follows from the fact that any asymptotic measurement will 
always be governed by projectors on some out-states, i.e. states with a definite particle content 
defined on J + , to probe the physics. Therefore, these measurements will always be of the 
in — out type since the Heisenberg state of the field is specified (prepared) before the collapse. 
It is this two-states formalism giving rise to non-diagonal matrix elements [5] (exactly like in 
a 5-matrix formulation [12]) which is at the origin of the difference: The metric fluctuations 
cannot coherently affect configurations specified in the 'kef on J~ and in the 'bra' on J7" , 
hence the coherence is lost. On the contrary, measurements performed by infalling observers 

17 One might wonder if the effects we are describing are not induced by the choice of working at fixed u or at fixed 
v, r. To avoid misunderstanding, we recall that Green functions have no physical meaning per se, rather they are 
elements which appear in integrals describing transition amplitudes (for a discussion of this point in a quantum 
gravitational context see Section 2 in [50]). Having made this remark, one verifies that u is a physically meaningful 
coordinate on J + by noticing that du\ r — At where At is the proper time of a particle detector on J + . Thus 
when additional quantum mechanical systems are coupled to the radiation field, the matrix elements governing 
transition amplitudes will have, in their integrand, phase factors behaving like e _!A11 in any coordinate system. 
Similarly, upon questioning what an infalling observer might see when crossing the horizon, the dependence in 
v,r is meaningful since Ar\ v oc — dr where dr is his proper time. 

18 This clearly illustrates that the physics seen by infalling observers completely differs from that reconstructed 
from observers at large distance from the hole. This is similar to what was advocated in [TJ|. However, the fact 
that the FT transformed on J + is governed by eq. (|43p indicates that the near horizon physics is unaccessible 
and therefore lost to remote observers. Thus it seems that these two descriptions cannot not obey the 'comple- 
mentarity principle' [TJ]. Indeed, it appears from our analysis that in the interacting FF vacuum early out going 
configurations are entangled to infalling ones. Let us recall that by complementary to each other, it was meant 
that the two descriptions are both complete, like the position and momentum representations of the same vector 
state in quantum mechanics. 
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only probe the near horizon behaviour of the Green function. Hence the coherence is maintained 
for them. 



6 Conclusions 

We have studied the effects induced by the gravitational interactions governed by eq. (|16p . Even 
though we worked out only the lowest order in G (cja oc G) we believe that our main result is 
robust We see no reason for higher order terms to suppress the entanglement of 0_ with 
so as to give a\ = thereby recovering trans-Planckian correlations. Indeed higher order 
modifications to eq. (|37p and eq. (|38p should be of the type {Guj 2 ) n and therefore will not 
affect the low frequency behaviour of eq. (|38p thereby leaving the effective spread d\ essentially 
untouched. Moreover, when considering the effects of higher angular momentum modes, as 
indicating by [15], a should be larger than our estimate based on s-modes because the effects of 
higher angular momenta should add incoherently. 

In brief, given that gravitational interactions grow without bound near the horizon, we claim 
that the entanglement of </>_ with <j) + is unavoidable and universal. (By universal we mean that a 
similar entanglement would be also found when considering the coupling of <p to other quantum 
fields such as massive ones.) The entanglement will then prevent the unbounded growth of 
frequencies encountered in the free field theory and will be accompanied by the reorganization 
of the description of vacuum. That is, when r — ► 2M, the usual states of the free field theory 
which give rise to the notion of on-shell asymptotic particles provide bad approximations of the 
true eigenstates (albeit still characterized by the Killing energy A since the situation is stationary) 
of the interacting the ory . It is the growing discrepancy which leads to the 'dissipation' of the 



We also claim that the radiative corrections encoding dissipation should be finite because 
only low frequency to ~ k infalling configurations contribute to them. Indeed, in the UV regime 
for both infalling and outgoing configurations, the expressions coincide with those evaluated in 
the tangent plane and are Minkowskian in character. Hence the high frequency regime cannot 
contribute to the effects which break Lorentz invariance. (This still needs to be confirmed by 
an explicit calculation.) 

We would like to conclude this work by several remarks on related aspects of Quantum 
Gravity and Black Hole Physics. 

First, we point out the similarity between the above effects induced by the metric fluctuations 
representing gravitational interactions in the vacuum and those attributed [51] to 'foam'. By 
foam we mean quantum gravitational configurations which radically affect the smoothness of 

19 Added comment : We can not exclude the possibility that, when properly renormalized, the 1-loop radiative 
corrections governed by fg6'|) will give a a = 0. In fact, the self-energy of an electron in a thermal bath of photons 
(as a thermal bath of gravitons) induces no dissipative effects at one loop, see [28\ and refs. therein. Nevertheless 
the naive reasoning that loop corrections in a thermal bath induce dissipative effects 'becomes' correct at two-loop 
level. These effects are indeed propotional to e 4 where e is the electric charge. Thus, in the case one finds <ta = 
at one loop, one still needs to confront the question whether this result be found at any loop level, to all order in G 
(for a deep reason, e.g. because we are dealing with the true interacting vacuum), or shall <ta 7^ at some higher 
loop level (as in a heat bath, because the integrands are no longer 2D Lorentz invariant), thereby validating the 
present conclusions. 

20 Motivated by the present work, we have recently constructed an interacting QFT model which is exactly solvable 
and which explicitly displays this phenomenon; '4-91, namely the growing discrepancy between the usual modes and 
the true eigenstates manisfests itself through the dissipation of the two-point function. 
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space-time at short distance. (They might arise from gravitational instantons [52] or stringy 
effects [S3].) In all cases, the replacement of free field propagation in a fixed background by 
the appropriate interacting model might lead to very similar (universal?) deviations when 
analyzing the departure from the free field description that all models possess at large distances. 
Therefore, these first deviations might be described by some effective mesoscopic theory of space- 
time properties which would essentially signal the existence of a minimal resolution length [54], 
the equivalent of our a, in the otherwise local field theory. 

Secondly, we conjecture that a (properly computed so as to include the contribution of higher 
angular momentum modes) should also be the length scale which enters in the entanglement 
description of the black hole entropy [55J. We recall that when using free field in a given space 
time, the entanglement entropy diverges due to the unbounded character of the reservoir of high 
energy modes. To get a finite entropy density per unit area, some cut-off should be introduced. 
We propose that the cut-off defining the black hole entropy should be the dynamically induced 
length-scale cr(iV), i.e., the length scale at which correlations between configurations on J + and 
the near horizon region get lost when N quantum fields contribute to the entropy and therefore to 
the near horizon gravitational interactions. We refer the interested reader to our recent work \27$ 
wherein we study the entanglement entropy in the presance of dispersion defined in a FF frame 
such as that of eq. {5]j. In that work, it is shown that one obtains an entanglement entropy which 
scales as the horizon area per Planck length square only when the UV cutoff ( which specifies at 
which scale measured in the FF frame the propagation ceases to be the usual Lorentz invariant 
one) scales as in our model, i.e. when the UV cutoff scales as 1/a cx Mp lanck / k. 
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7 Appendix: The large-iV limit 

We briefly mention several interesting features of the large N limit which illuminate the problems 
we addressed. 

The semi-classical description of quantum processes occurring in a curved background is 
based on the following equations 



In eq. (j49l) the field operator propagates in the classical geometry g = g* which is a solution 



G, 



8ttG(*|T^|*) 



(48) 
(49) 



U g( j> = 0. 
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of eq. (j35J) driven by the expectation value (^T^l^) evaluated in the state |^) using eq. (09]). 
In this sense, eq. ([151 |4*9"P is a self-consistent (Hartree) approximation. 

It is quite reasonable that this approximation correctly predicts the time evolution of certain 
quantities in certain circumstances, e.g., the rate of mass loss of a large black hole. However, the 
criteria which characterize the validity range of the predictions obtained from eq. (|48l49p are 
not known. An obstacle in finding these criteria is the identification of the "small parameter(s)" 
which control the deviations between the exact evolution and the semi-classical one 

A rather formal answer to these questions is provided by considering a large N limit, where N 
is the number of copies of the <f> field. The simplest way to understand why the above equations 
govern the large N limit is by considering the path integral approach of matrix elements. By 
duplicating N times <fi in eq. (Ilip and first integrating over them with a fixed metric h, the 
1-loop effective action for h contains an overall prefactor N when replacing Newton's constant 
G by G/N . Therefore, in a large N limit, the path integral over h can be evaluated by a saddle 
point approximation. The stationary phase condition then gives rise to eq. ([48 p where the 
expectation value of Tn V is that of one field. The spread around the saddle point scales like 

AT-l/2. 

In this approach, the validity of the semi-classical equations apparantly relies on a statistical 
argument, as mean quantities emerge in the thermodynamic limit. The weakness of this argu- 
ment is the absence of role played by the hierarchy of the length scales governing the processes 
under examination. This is unlike what is found when a Born-Oppenheimer treatment is applied 
to quantum gravity |58j. 

More interestingly the large N limit makes also predictions beyond the semi-classical equa- 
tions. For instance, in the limit N — > oo with GN fixed, the short distance behaviour of the 
graviton propagator is modified, see |43 [ 1171 HI] . In Minkowski vacuum, these modifications are 
of course Lorentz invariant. However, in a thermal bath or a curved background, the correction 
terms will no longer possess the Lorentz invariant form. Hence they can induce the effects we 
are seeking: a dynamically induced scale which breaks the (local) Lorentz invariance that the 
un-interacting theory possessed. Moreover, only low energies {i.e., energies comparable to the 
temperature) contribute to this new scale because in the UV limit all expressions tend to their 
Minkowski vacuum, and hence Lorentz-invariant, form. Hence they shouldn't be any additional 
UV divergences in the expressions giving rise to the new scale, as it is the case for the corrections 
to the self-energy of an electron immersed in a thermal bath. 

To further strengthen the relations with what we did in Section 6, it is instructive to see 
how the semi-classical treatment emerges from eq. (j35p viewed as generating perturbatively the 
connected Feynman diagrams when expanding e lSint in powers of Si n t- In this description, one 
finds that the Green function is a double sum of powers of N and G which possess the following 
properties. 

• The power of N is equal or inferior to that of G. 

• The semi-classical treatment corresponds to the leading series: the set of graphs weighted 
by (GN) n . All graphs are one-particle reducible and are governed by the one-point function 
(Tu U ). Upon summing this series, one identically recovers the Green function evaluated in 
the 'mean' geometry the solution of eq. (j48[) . 

• Having summed up the leading series in (GN) n , our treatment corresponds to the next 
series: the set of graphs weighted by (G 2 N) n . All graphs are two-particle reducible and 
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are governed by the connected two-point function (T^Tap) c . Upon summing this new 
series, one obtains the Green function evaluated in the stochastic ensemble governed by 
eq. (f38j) . This second series should also be related to the use of the above mentioned 
large- N modified graviton propagator in the place of the unperturbed one. 

In brief, N is a parameter which organizes the double sum of graphs into a sum of non- 
perturbative series. When the former series have been summed up, the m-th series contains 
all powers of (G- m iV), and is governed by the 77?*-th correlation function of T^. 

The physical question raised by these results is the following: given the dimensionality of 
G = lpi anck , can one infer that high orders in m (the relative power of G with respect to that 
of iV) become relevant only for high (Planckian) energies ? We conjecture that this is the case: 
the sorting out of graphs in terms of m is effectively an expansion in the energy of the processes 
involved in the matrix element under consideration. This is what seems to emerge from our 
analysis. As indicated by eqs. (|431 144p . the semiclassical description of the correlations breaks 
down when r — 2M ~ a, i.e., when the FF energy of a mode f2 ~ p r = \/(r/2M — 1), see eq. 
(|42p . reaches the new scale l/crffil We thus find, as in |58j . that the validity of the semi-classical 
equations relies on a hierarchy of energy scales. Indeed, for a large black hole, a <C M even 
when N = 1. Thus N is not necessary to justify the validity of the semi-classical description. 
It is rather a useful parameter which helps sorting out the different contributions in radiative 
corrections. In particular, it allows to discard the self-interacting (non-planar) graphs which are 
more difficult to evaluate (since they cannot be expressed in terms of metric fluctuations in the 
present case) because their power of N is smaller those weighted by (G m N) n . 
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